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Abstract. We report the results of an observation with the Rossi X-ray Timing Explorer of the black hole candidate 
GX 339-4 during its 2002/2003 outburst. This observation took place during a spectral transition from the hard 
to the soft state. A strong (6% rms) transient quasi-periodic oscillation (QPO) appears suddenly in the power 
density spectrum during this observation. The QPO centroid is ~6 Hz, but it varies significantly between 5 and 

7 Hz with a characteristic time scale of ~10 seconds, correlated with the 2-30 keV count rate. The appearance of 
the QPO is related to spectral hardening of the flux, due to a change in the relative contribution of the soft and 
hard spectral components. We compare this peculiar behavior with results from other systems that show similar 
low frequency QPO peaks, and discuss the results in terms of possible theoretical models for QPO production. 
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1. Introduction 

Variability studies of black-hole candidates (BHCs) have 
revealed the presence of many types of quasi-periodic 
oscillations (QPOs ) in the X-ray flux of these systems 
(van der Klislll995t iRemillard et all l2002al) . While most 
BHC have shown QPOs between 0.01 and 30 Hz, in sev- 
eral of them also QPOs above 40 Hz have been detected. 
In particular, QPOs between 1 and 10 Hz are very com- 
mon in these systems and are p robably linked to sim ilar 
features in neutron star systems ijBelloni et aljE00 2b1 . As 
both the low and high frequency QPOs are thought to 
arise in the accretion flow close to the black hole, they are 
a potentially important tool in determining the properties 
of these compact objects and the accretion flow. 

The presence and the properties of QPOs in BHCs 
seems to be related to the spectral characteristics of the 
source. In this respect, transient BHCs are particularly im- 
portant for the study of QPOs, since these systems show 
outbursts during which the mass accretion rate changes by 
several orders of magnitude, usually resulting in several 
transitions between different accretion modes or states. 
These states are characterized by distinct spectral and 
temporal properties. QPOs are mostly found in states in 



which the hard spectral component contributes signifi- 
cantly to the spectrum and they are usually stronger at 
higher energies. Although a considerable amount of data 
showing QPOs has been collected over the years, many is- 
sues such as the nature of the frequency changes, appear- 
ance/disappearance of the QPOs, relation to less coherent 
variability, etc., have remained unsolved. Several theoreti- 
cal models exist for these oscillatio ns, but there is no con- 
sensus as to their physical nature lAlpar fc Shahamll 19851: 
Strohmaver et al.lll996t iMiller et al.lll998t IStella fe Vietril 
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1998) . 

GX 339-4, for whic h the recently est ablished mass 
function (5.8 ±0.5 M^. lHvnes et al.l (2003^1 indicates the 
presence of a black-hole primary, is usually referred to as 
a persistent BHC. However, its long term behavior is more 
similar to that of transients and it is one of the few sources 
that has shown all black hole spectral /timin g states 
(see iTanaka fc Lewinl fl995t iHoman et all l200ll iNowakl 
2002J for a detailed description of black hole states). 
While GX 339-4 is usually observed in the low/hard 
state (LS), several state-transitions have been reported 
llMaeiima et al.lll984tlMivamoto et"aIlll99ltlBelloni et alJ 

1999) . Strong QPOs were observed in the power den- 
sity spectra of GX 3 39-4 during its state transitions. 
iMivamoto et all l)l99l(l reported 1-15 Hz QPO together 
with strong variable band-limited noise in the VHS tran- 
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sition observed with Ginga. A complete analysis of simi- 
lar osciJlMionsj^ GS 1124-68 was performed 
bv lTakizawa et al who showed that there are two 
different sets of QPOs, with different dependence on the 
source count rate. During these observations, fast count 
rate variations were observed, coin ciding with changes in 
the spectral and timin g properties ^Miyamoto et alJfl99lt 
iTakizawa et al.lll997|) . So far, no high frequency QPOs 
have been detected in GX 339-4. 

After its 1998/1999 outburst, GX 339-4 returned 
to the LS, and it became undetectable for the All- 
sky Monitor (ASM) onboard the Rossi X-Ray Timing 
Explorer (RXTE) in May 1999, entering its 'off' state, 
obser ved before only once with EXOSAT ijllovaiskv et al.1 
1986). The source was detected with BeppoSAX on 
1999 August 13, and its properties were consistent 
with th e off state being a low-luminosity extension o f 
the LS l|Mendez fc van der Klislll997t iKong et all l2000l) . 
However, a later observation with BeppoSAX on 2000 
September 10 revealed a flux level l ower by a factor o f 
~4 than what was detected in 1999 l|Corbel et alJl200.^ . 
GX 339-4 has remained undetectable by the ASM un- 
til 2002. In March 2002, the source became bright again 
llSmith et al.i r20021 and we started an RXTE campaign 
to follow its evolution. Here we present the results of the 
timing analysis of a single observation from this campaign, 
in which a QPO with exceptional behavior was observed. 
The full spectral/timing analysis from all observations of 
our campaign will be presen ted elsewhere. Prelim inary re- 
sults have been reported by Bclloni et alJ l(2002a^ . 
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Fig. 1. The 2.9-25.1 keV count rate (top) and color (bot- 
tom) evolution of GX 339-4 during its 2002 /2003 outburst , 
as observed with RXTE (PCU2 only). The color is defined 
as the ratio of counts in the 16.4-19.8 keV band over counts 
in the 3.7-6.5 keV band. Day corresponds to April 3 2002 
(MJD 52367), the vertical dotted line marks the time of 
the observation analyzed here (May 17 2002/MJD 52411, 
see Fig. 2 for blow up), and the gray area indicates the 
transition interval. The time resolution of the light curve 
is 16s, while in the color curve each point represents a 
single observation. Data are corrected for background but 
not for dead time. 



2. Data analysis 

Figure □ shows the 2002/2003 outburst of GX 339-4. The 
light and color curves are obtained from the data of the 
Proportional Counter Array (PCA) onboard RXTE. Only 
counts from the most reliable of the five units of the 
PCA (PCU2) were used. From the figure it is evident 
that the outburst had a complex evolution, which can be 
roughly divided in three parts. In the first part (day 0- 
34), the count rate increased with time, showing a rather 
hard spectrum gradually softening with time. The second 
part (day 34-57, depicted by the gray strip) showed a 
drop in count rate accompanied by considerable spectral 
softening. During the third part (day 57-292) the count 
rate returned almost to the level of the first peak, after 
which it showed an overall, but sometimes irregular, de- 
cay. These changes occurred while the spectrum remained 
much softer than during the first part of the outburst. 
From the light curve, it is also clear that variability on a 
time scale of tens of seconds is much stronger in the spec- 
trally hard part of the outburst. Based on the timing and 
color data, we conclude that the source evolved from the 
LS, through the very high state (VHS) o r intermediate 
state (IS), to the HS l)Belloni et al.ll2002a|) - A dotted line 
in Fig. ^ marks the observation discussed in this work. 

This observation took place during the transitional 
phase of the outburst and, although not outstanding in 



the light curve, showed peculiar timing features, differ- 
ent from those observed before and after. It corresponds 
to 2002 May 17/MJD 52411 (two RXTE orbits, 14:24- 
14:58 UT and 15:34-16:33 UT), with a total PCA exposure 
time of 5.58 ks. In Fig. [2 we show a 16-s time resolution 
light curve and the color evolution during this observation. 
Below, we present the results of the timing and spectral 
analysis from these data. 

2.1. Power Density Spectra 

For the timing analysis, we used Single bit data from 
all active PCUs, with a time resolution of 1.25xl0 _8 s 
and covering the PCA channel range 0-35, corresponding 
roughly to 2-16 keV. For each of the two orbits, we sub- 
divided the data into stretches 128 seconds long and pro- 
duced a power density spectrum (PDS) for each of them. 
The resulting PDS were then averaged, a logarithmic re- 
binning was applied, th e contribution due to Poissonia n 
statistic was subtracted ljZhandll995l IZhang et alJll995|) . 
and finally they were conv erted to fractional squared rms 
llBelloni fc Hasingedll990T) . The resulting power density 
spectra for each orbit are shown in Fig. |3| The two PDS 
are markedly different. In the first orbit, the power density 
spectrum shows weak power-law noise with a very broad 
bump between 6 and 10 Hz. In the second orbit, the power- 
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law noise is enhanced and a strong QPO at ^6 Hz appears, 
together with its second harmonic. We fitted both power 
density spec tra with a model consisting of Lorentzia n 
components l|Mivamoto et alJll99lt iBelloni et al.ll2002bj) . 
During the first orbit, the broad QPO (centroid frequency 
7.3±0.2 Hz and full-width at half-maximum 4.4±0.4 Hz) 
and the continuum noise (consisting of two Lorentzian 
components) have a fractional rms of 1.73±0.09% and 
1.92±0.11% respectively. A Lorentzian function proved to 
be a poor approximation for the strong QPO in the sec- 
ond orbit. This effect has also been observed in a QPO 
at a similar frequency in XTE J 1550-564 (Homan et al. 
2001). A good fit to the QPO peak was obtained with a 
model consisting of a Gaussian plus a broad Lorentzian to 
account for the extended wings. The centroid frequency of 
the Gaussian -(-Lorentzian complex is 5.895±0.004 Hz (the 
values for the two components are constrained to be the 
same). The width of the Gaussian is 0.351±0.006 Hz, and 
the total (Gaussian plus Lorentzian) fractional rms ampli- 
tude is 6.68±0.07%. The continuum noise, fitted with two 
broad Lorentzian components (as for the first orbit), has a 
total fractional rms of 5.06±0.08%. We repeated the anal- 
ysis with data at higher energies (PCA channels 36-51, 
corresponding to 17-24 keV). No signal is detected dur- 
ing the first orbit, although with not very stringent upper 
limits. During the second orbit, the QPO fractional rms 
is 17.81±0.36%. 
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Fig. 3. Power Density Spectra from the two RXTE orbits 
(top: first orbit; bottom: second orbit). The solid lines 
correspond to the best-fit models described in the text. 
For reference, the dashed line in the lower panel shows 
the best fit of the first orbit. 



centroid frequency, we produced a spectrogram, defined 
(in its continuous version) as 
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Fig. 2. A blow up of the light and color curve for MJD 
52411. Time resolution is now 16s in both curves. Typical 
error bars are shown. See caption of Fig. 1 for additional 
details. 



2.2. Time-frequency analysis 

In order to test whether the inadequacy of the Lorentzian 
model for the QPO is due to intrinsic variations of its 



Here, x(t) is our signal (PCA channel range 0-35, with 
a resolution of 128 _1 s) and h(t) is a window function 
(in our case, a W elch window of length 4 seconds, see 
IPress et al.l ll992). The time step t, which sets the time 
resolution of the resulting time- frequency image S x (t,u) 
was chosen to be 1 second. In order to suppress the low- 
frequency noise, before applying Eqn. 1 we detrended our 
signal by subtracting from it a spline fit to its 1-s rebinned 
version. A section of the spectrogram, showing 320 seconds 
(starting from 50 seconds into the second RXTE orbit) in 
the frequency range 4-8 Hz is shown in Fig. 0] together 
with the corresponding light curve (in 1-second bins). Four 
effects can be observed: 

— The QPO is not present in the first 120 seconds of the 
observation, (we can set a 95% confidence upper limit 
of 2.1% rms by fitting the PDS of the first 120 s with 
a QPO with centroid and width fixed to the values 
reported above), barely visible in the next 20 seconds, 
and appears clearly only after that (to remain for the 
rest of the RXTE orbit) . 

— The light curve shows a net increase in count rate in 
correspondence with the onset of the oscillation, when 
also the enhanced low-frequency noise, visible in Fig. 
13 turns on. 

— The centroid frequency of the QPO shows significant 
variations in the range 5.5-7.0 Hz. 

— The centroid frequency shows a positive correlation 
with the count rate (shown in the top panel). 
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Fig. 4. Top panel: PCA light curve (1 second binning) 50- 
370 seconds into the second RXTE orbit. Bottom panel: 
corresponding spectrogram in the 4-8 Hz range; darker 
regions correspond to higher power. 



On the basis of third point we conclude that the shape 
of the QPO, in particular its broad wings, are indeed the 
result of the variations in the centroid frequency. 

To examine in more detail the QPO-count rate rela- 
tion, we fitted each power density spectrum in the spec- 
trogram (in the 4-8 Hz range) with a model consisting 
of a Lorentzian peak and a power law for the local con- 
tinuum. Visual inspection of the results showed that the 
model was able to recover accurately the local centroid fre- 
quency of the QPO. On average, the FWHM of the best 
fit Lorentzian was found to be ~0.2 Hz. A correlation be- 
tween QPO frequency and count rate over the entire orbit 
is clear, although it shows a large scatter. 

Since the QPO centroid frequency varies in an irregu- 
lar way, we produced a PDS of the QPO-frequency time 
series. In order to avoid problems due to the fact that adja- 
cent Welch windows overlap with each other, introducing 
a net filtering effect, we repeated the whole analysis with 
a window length of 2 seconds and non-overlapping win- 
dows. The resulting power density spectrum is shown in 
Fig.|SJ A break in the power density spectrum is evident. A 
fit with a zero-centered Lorentzian yields a characteristic 
frequency (HWHM) of 0.09 Hz. Fig. [5] shows the corre- 
lation between QPO frequency and total count rate from 
these 2-second windowed data. One can see that there is 
indeed a correlation between the two quantities, but it 
is not very strict, with a linear correlation coefficient of 
~0.5. Interestingly, the additional noise component below 
0.1 Hz in the second orbit, reflects similar time scales as 
the variations of the QPO frequency. 

Although variability properties, such as QPO frequen- 
cies, generally correlate better with spectral hardness than 
with count rate, the statistical quality of our data pre- 
vented us from performing a meaningful study with a color 
instead of count rate. 



Table 1. Best fit parameters from the spectral fitting of the 
data from the two RXTE orbits. Errors correspond to a la 
confidence range. The total unabsorbed 2-100 keV luminosity 
was c alculated assuming a distance of 4±1 kpc iZdziarski et alJ 
1998). 



Orbit 1 Orbit 2 



PL Index 


2.44L0.02 


2.44±0.10 


E c (keV) 


1901JS 


11 0+ 70 


T in (keV) 


0.89±0.01 


0.92±0.01 


E ime (keV) 


6.86L0.11 


6.7±0.2 


F To t a 


(2.06±0.06) x 10 -8 


(2.25L0.17) x 10" 8 


"FT a 
-T disk 


(1.30±0.04) x 10~ 8 


(1.18L0.02) x 10" 8 


F a 
A pow 


(7.71±0.35) x 10~ 9 


(1.02±0.03) x 10~ 8 


L>tot (erg s" 1 ) 


(3.73L1.86) x 10 37 


(4.07±2.05) x 10 37 



a unabsorbed 2-100 keV flux (erg cm 2 s 1 ) 



2.3. Energy spectra 

We performed a spectral analysis of the two individual 
orbits separately, ignoring the first 120 s of the second 
one. Spectra were created from PCU2 data (2.9-25.2 keV) 
and HEXTE cluster A data (18.3-150 keV), using the 
latest version of FTOOLS (5.2). The spectra were cor- 
rected for background and dead time effects, and fitted 
with XSPEC 11.2 using a model consisting of a cutoff 
power law (cutoff), a disk black body (diskbb), an ac- 
cretion disk line (laor) and a smeared edge (smedge). The 
interstellar absorption Nh was fi xed to a value of 5 10 21 
atoms/cm 2 l)llovaiskv et aljfl98filh For both orbits, good 
fits were obtained. The fit parameters were consistent be- 
tween the two orbits (see Table 1). The only differences 
were found in the fluxes of the disk black body and of the 
cutoff power law: while the first decreased by 9%, the sec- 
ond increased by more than 30%, resulting in an overall 
hardening of the spectrum. The calculated fluxes in the 
2 - 100 keV range for the two RXTE orbits are reported 
in Tabic 1. 

3. Discussion 

We have discovered a strong transient 6 Hz QPO in 
GX 339-4 with some unusual properties: a sharp turn-on 
in time and large frequency variability on a typical time 
scale of ^10 seconds. The QPO was found in an observa- 
tion that was performed during a transition from a hard 
to a soft state, during which the spectral and variability 
properties were those typical for a Very High State. It 
is interesting to note that, while the observation occurred 
during a transition from a hard to a soft spectral state, the 
appearance of the QPO in our observation is associated 
with the hardening of the energy spectrum. This indicates 
that such a hard to soft transition is not a smooth process 
throughout the outburst. 

The sharp 6 Hz QPO in the second orbit turned on 
within ^20 seconds, similar to the timescale of the sharp 
increase observed in the light curve. While a broad QPO 
was present during the first part of our observation, it is 
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Fig. 5. Power Density Spectrum of the QPO centroid fre- 
quency (see text). Power is in units of Hz 2 
. The line is a best fit with a zero-centered Lorentzian model. 



not clear whether the 6 Hz QPO evolved from this broad 
feature; the former is much weaker, less coherent, and 
peaks at a different frequency. We note that in case the 
broad feature is still present in the second part of the ob- 
servation, it might be partly responsible for the asymmet- 
ric wings of the 6 Hz QPO (see Fig. EJ- Such 'shoulders' 
have been reported before for similar 6 H z QPOs in other 
BHCs . In t he terminology in t roduce d bv lWiinands et alJ 
<jl999h and iRemillard et aD l|2002bh for QPOs in XTE 
J1550-564, the peak in the first orbit could be interpreted 
as a type A QPO, and the ones in the second orbit as 
type B; the observed spectral hardening is consistent with 
this interpretation. Besides the changes around 6-7 Hz, 
also the continuum at lower frequencies shows consider- 
able changes. As noted before, the increase in power below 

0. 1 Hz might reflect the count rate variations that corre- 
late with the changes in the QPO frequency that have a 
similar time scale (see Fig.^Jl. 

The observed FWHM of ^0.2 Hz in the dynamical 
power spectrum is consistent with being solely due to 
smearing by frequency changes. This sets a lower limit 
to the Q value of the QPO of ~30, indicating a signal life- 
time of at least 5 s. The fact that we see uninterrupted fre- 
quency changes of 10 seconds and more suggests a longer 
lifetime. If we associate the observed frequency with a cer- 
tain radius in the disk, it is unlikely that the QPO is 
caused by a single blob of matter, as the blob would have 
to move inward and outward. It is possible that what we 
observe is the result of different blobs that appear and 
live at slightly different radii, which vary on a similar time 
scale. However, Fig. 01 and Fig. 5 show that the QPO fre- 
quency variations are not consistent with a white noise, 

1. e. there is a correlation between the QPO frequencies 
at different times. This would mean that each new blob 
would "know" about the previo us ones. In a model like 
that of IPsaltis fc Normanl l|2003h . where a particular ac- 
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Fig. 6. QPO frequency vs. PCA count rate for the 2-s 
binned data (see text). Typical error bars are shown. 



cretion disk radius acts as a band-pass filter, the observed 
variations would translate into corresponding variations 
of this filtering radius, resulting in a natural correlation 
between different frequencies at different times. The ob- 
served count-rate/frequency correlation would then sug- 
gest a relation between this filtering radius and the X-ray 
flux. 

Given its peculiarities, one might wonder whether this 
QPO detection is unique. A similar sharp QPO peak, 
ne eding a Gaussian mo del for a satisfactory fit, was found 
by iHoman et alJ l|200lf) in the black-hole t ransient XTE 
J1550 -564. This QPO (defined as type B bv lHoman et alJ 
l200lh was also detected in a few observations with a cen- 
troid frequency around 6 Hz. A sub- harmonic and two 
higher harmonics were also detected. Interestingly, a sim- 
ilar sharp QPO was observed in the Ginga data of the 
trans i ent GS 1124-68 ilMivamoto et alll994UBelk>ni et alJ 



Il997j) . iTakizawa et alJ l)1997|) classify the QPOs observed 
in GS 1124-68 into two separate classes: the sharp QPO 
peak was only observed in the range 5-7 Hz , simil ar to 
what we detected in our data. IHoman et al.l l(200l|) also 
report the sudden appearance of a QPO in XTE J 1550- 
564 associated to a count rate increase. However, in this 
case the X-ray colors indicate a small and gradual soften- 
ing of the spectrum. A sharp QPO around 6 Hz, with 
harmonic and sub- harmonic com ponents and a similar 
low-frequency noise, was found bv lCui et al.l 1 2000l) in the 
transient XTE J1859+226. A narrow peak was present 
in three of the four observations analyzed, with a cen- 
troid frequency around 5.96 Hz, with properties similar to 
those seen in our data. Another narrow QPO peak was 
foun d in the bright transi ent GRS 1739-27, although at 5 
Hz l|Wiinands et al.ll200l|) . The QPO, with its first over- 
tone, was observed in the second part of the data and not 
in the first part, where onl y a weak noise compon ent at 
low frequencies was noticed. IWiinands et a found 
that the frequency of the QPO increased together with the 
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count rate, with an even better correlation with the hard 
color. All these QPOs have a centroid frequency of 5-7 Hz. 
It is not clear, however, why the same frequency should 
play the same role in different systems. In all these cases, 
the QPO is associated to a VHS-like energy spectrum, 
with the presence of both a hard and a soft component, 
has a fractional rms of a few percent and a FWHM sim- 
ilar to that shown here. In order not to vary too much 
over systems with black holes of different masses, the de- 
pendence of its frequency with mass of the central object 
must be weak. 

The QPO observed here shows, moreover, characteris- 
tics in common with both the HBO and the NBO of NS 
LMXB (see van der Klis 1995), letting open the question 
if a stronger connection can be established either with one 
or the other feature. HBOs present LFN (Low Frequency 
Noise), with rms of 1-8 %, the oscillation has a frequency 
that correlates with count rate and an rms of 2-8 %, the 
spectrum is hard. NBOs show VLFN (Very Low Frequency 
Noise), with a rms of 0.6-4 %, no LFN, the frequency of 
the oscillation varies between 5 and 7 Hz, and its rms is 
1-3 %, while the energy spectrum is hard. Nothing can be 
said about correlation between count rate and frequency 
over short time scale. 

The rapid count-rate increase observed in our ob- 
servation, corresponding to a hardening of the flux, 
shows sim ilarities to the "dips" and "flip-flops" re- 
po rted bv iMivamoto et all l)l99ll) from GX 339-4 and 
bv lTakizawa et alJ l|l997l) from GS 1124-68, although the 
time scale of our transition is slower. In both cases, 
GX 339-4 and GS 1124-68, a QPO peak was observed only 
in the high-flux intervals, similar to our data, although 
the continuum noise was stronger at low flux, contrary to 
what we observe here. The corresponding spectral varia- 
tions were similar to ours: the flux increase was caused by a 
30% increase of the hard component, while the disk com- 
ponent did not vary. The spectra l parameters remained 
unchanged l|Mivamoto et alJll99l|) . Notice that also in GS 
1124-68, the centroid frequency of the sharp QPO is pq s- 
itively correlated with count rate i Takizawa et alJ[l997h . 

The spectral hardening observed in connection to the 
appearance of the QPO suggests that the oscillation is as- 
sociated to the hard spectral component, as it is usually 
observed for low-frequency QPOs. However, the presence 
of a soft component seems to be necessary for the obser- 
vation of such a QPO, as one can see also in the cases of 
XTE J1550-564, GS 1124-68 and XTE 1859+226. A QPO 
with these characteristics (narrow, around 6 Hz and with 
a weak noise component) seems to be always associated to 
a sort of "intermediate" spectral state, when both spec tral 
components are present (see e.g. iRutledee et aklll999|) . 

The results showed here pose constraints to theoretical 
models for the production of the QPO. A successful model 
should be able to reproduce a sharp appearance of a QPO 
like the one we observe, but also the rapid variability of 
the characteristic frequency of the oscillation (see Fig. EJl 
should be explained. 
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